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Abstract 

Spectral properties for the two-dimensional quantum S = 1 XY model were 
investigated with the exact diagonalization method. In the symmetry-broken 
phase, there appear the massive Higgs and massless Goldstone excitations, 
which correspond to the longitudinal and transverse modes of the sponta¬ 
neous magnetic moment, respectively. The former excitation branch is em¬ 
bedded in the continuum of the latter, and little attention has been paid to 
the details, particularly, in proximity to the critical point. The finite-size- 
scaling behavior is improved by extending the interaction parameters. An 
analysis of the critical amplitude ratio for these mass gaps is made. 
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1. Introduction 

In the symmetry-broken phase, the 0(2)-symmetric system, such as the 
XY model, exhibits a massless Goldstone excitation, which corresponds to 
the transverse modulation of the magnetic moment. On the one hand, the 
longitudinal mode, namely, the Higgs excitation, is massive, embedded in the 
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continuum of the former; see Ref. [lj for a review. The 0(2)- [equivalently, 
U(l)-] symmetric system is ubiquitous in nature, and such a characteristic 


spectrum has been observed for a variety of substances 


, la, m la, m, la, Is, 


3, llO, [nj. The perturbation he 


0(2) (axial) symmetry 




d (experimental probe) should retain the 


141, Il5| : otherwise, the contri 


the Goldstone excitations smears out the Higgs-mode branch 


16 


jut ion from 


(For instance, the chemical-potential modulation for the bosonic system does 
not conflict with the symmetry.) 

Recent studies 2Q, 2l| shed light on a universal character of the spec¬ 
trum in proximity to the phase transition, especially, in (2 + 1) dimensions; 
it would be intriguing that the spectral property is also under the reign of 
universality. In (3 + 1) dimensions, the criticality is described simply by the 
Ginzburg-Landau theory (Gaussian fixed point). On the contrary, in (2 + 1) 
dimensions, the spectral property is non-perturbative by nature. In partic¬ 
ular, a universal amplitude ratio for the mass gaps [see Eq. (JT|) mentioned 
afterward] is arousing much attention recently. 

In this paper, we investigate the two-dimensional quantum S = 1-spin 
XY model ([2]) with the exact diagonalization method. The method enables 
us to calculate the low-lying level indexed by quantum numbers. In order 
to suppress corrections to scaling, we incorporate various types of interac¬ 
tion parameters in addition to the ordinary nearest-neighbor ferromagnetic 
interaction Jnn- Thereby, we investigate the universality for the critical 
amplitude ratio 


m H /A = 


m H (Jj 


NN 


( 1 ) 
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with the Higgs mass rn h , the Goldstone mass me, and the reflected gap 
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A = jtiq(2J^n — Jnn~) with respect to the critical point technical details 


and underlying physics are explained in t 
has been estimated as mu /A = 2.1(3) 


ne nex 
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;■ section. The amplitude ratio 


23 1 and 3.3(8) 24] by means of 


the (quantum) Monte Carlo method. According to the recent_claborated 


25], 2.2 26J], 


renormalization-group analyses, the ratio was estimated as 2.4 
and 1.67 [2; 

To be specific, we present the Hamiltonian for the S = 1 XY model 


nj _ t \ ' (qx nx I qy qy\ j \ ' (qx ox < cv qy\ 

n — —JnN / ‘-’j + ^j ) — J NNN 2_ ^ Cl ^j + *- 5 i ) 

(ij) «h» 

+On Y, (Sf + S-+SI + Sff + dY(S*?- (2) 

[ijkl] i 

Here, the quantum S = 1-spin operator S* is placed at each square-lattice 
point i. The summations, p ))i an d Xqp'wp run over all possible 

nearest-neighbor, next-nearest-neighbor, and plaquette spins, respectively. 
The parameters (Jnn, Jnnn, Da) are the corresponding coupling constants. 
The parameter D denotes the single-ion anisotropy. We survey the coupling- 
constant subspace 


(Jnn, Jnnn , Da, D ) = ( Jnn, JnnnJnn/ Jnn , Da, D*), (3) 

parameterized by Jnn- At Jnn = Jnn> th e system undergoes a phase tran¬ 
sition; a schematic phase diagram is presented in Fig. [0 Here, the critical 
point 


(J* NN , J* NNN i D 'n, D *) = (0.158242810160, 0.058561393564, 0.10035104389, 0.957), 

(4) 


was adjusted 


28] to an IR fixed point with almost eliminated irrelevant 


in¬ 


teractions; that is, the coupling constants (Jnni Jnnni Da) were determined 
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through an approximative real-space renormalization group, and the remain¬ 
ing one D* was finely tuned via the conventional finite-size-scaling analysis. 
As shown in Eq. (l2lh the S = 1-spin model allows us to incorporate various 
interactions such as the single-ion anisotropy, with which one is able to re¬ 
alize the XY -paramagnetic phase transition. In this sense, the extention of 
the magnetic moment to S' = 1 is essential in our study. 

The rest of this paper is organized as follows. In Sec. [2j we present the 
simulation results. Technical details are explained as well. In Sec. El we 
address the summary and discussions. 


2. Numerical results 

In this section, we present the simulation results. To begin with, we 
explain the simulation technique. 


2.1. Simulation algorithm 

In this section, we explain the simulation algorithm. As mentioned in 
Introduction, the XY model (J2]) was simulated with the exact diagonalization 


method. We implemented the screw-boundary condition 


291] in order to 


treat a variety of system sizes N = 10,12,..., 22 (AT: number of constituent 
spins) systematically; note that conventionally, the system size N is restricted 
within N = 9,16,.... We adopt the algorithm presented in Sec. II of Ref. 
. The linear dimension L is given by L = \fN ; note that the N spins 


constitute a rectangular cluster. 

Thereby, we evaluated the mass gaps m# and mg via the following 
scheme. The exact diagonalization method yields the low-lying energy levels 
E 0 < Ei < E 2 < ... explicitly. Each level E t (k, S^ ot ) is specified by the 
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momentum k and the perpendicular magnetic moment S ^ ot ; in practice, the 
numerical diagonalization was performed within the subspace (k,S^ ot ). The 
Higgs- and Goldstone-mass gaps are characterized by 


m H = Ei(0,0) - E 0 (Q, 0), 

(5) 

m G = E 0 ( 0,1) - £0(0,0), 

(6) 


respectively. The reflected gap A with respect to the critical point J^ N is 
calculated by A = m G (2Jf [N ~ Jnn)- The gap m G (A) becomes massive in 
the paramagnetic (XY) phase Jnn < Jnn ( Jnn > Jnn )> an d hence, the 
ratio millX makes sense in the XY phase. The gap A > 0 is interpreted as 
the insulator gap through regarding the ladder operators Sf as the bosonic 
creation-annihilation operators. 

2.2. Scaling analyses of m HG 

In this section, we investigate the scaling behaviors for the mass gaps m# 
(ED and m G (ED- 

In Fig. [21 we present the scaled Higgs gap Lm# for the nearest-neighbor 
ferromagnetic interaction Jnn and various system sizes N = 10,12,..., 22 
(L = y/N). The data merge around Jnn ~ 0.15, indicating an onset of 
criticality; note that the scaled energy gap should be scale-invariant at the 
critical point. The location of the critical point is consistent with Eq. (1T|) . 
The Higgs gap mu appears to open in both the paramagnetic (Jnn < Jnn) 
and XY (Jnn > Jnn ) phases; the latter case is our main concern, as 
mentioned in Introduction. 
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In Fig. El we present the scaled Goldstone gap Line for various J^n 
and N = 10,12,..., 22. The data indicate an onset of criticality around 
Jnn ~ 0.15. The Goldstone gap closes in the XY phase, while it opens 
in the paramagnetic phase. The latter gap is interpreted as the (bosonic) 
insulator gap A, setting a fundamental energy scale in this domain; actually, 
the transition is interpreted as the superfluid-insulator transition 2l|. Hence, 
the ratio m#/A, Eq. (JT]), makes sense in the XY phase, and the criticality 
is investigated in the next section. 

In Fig. O we present the scaling plot, (Jnn — Jn^^^-Luih, for var¬ 
ious system sizes N = 10,12,..., 22. Here, the scaling parameters are set 


to J* NN = 0.158242810160 [Eq. g}] and v = 0.6717 (30!, |3lj. The data 
appear to collapse into a scaling curve satisfactorily. Similarly, in Fig. El 
we present the scaling plot, ( Jnn — Jnn)^ v ~L m G, for various system sizes 
N = 10,12,..., 22; the scaling parameters are the same as those of Fig. [4j 
We address a few remarks. First, the data in Figs. [I] and 0 collapse into 
the scaling curves satisfactorily. Such a feature indicates that corrections to 
scaling are almost negligible owing to the fine adjustment 28] of the coupling 
constants to Eq. g]). Because the tractable system size with the exact 


diagonalization method is severely restricted, it is significant to accelerate 
the convergence to the scaling limit. Second, the scaling parameters, J^ N 
and v, are taken from the literatures, Refs. 28] and 30], respectively. That 
is, there are no adjustable ad hoc parameters in the present scaling analyses. 
Last, as demonstrated in Figs. [2] and [3] both mass gaps rriuxj possess an 
identical scaling dimension. Hence, the amplitude ratio g]) makes sense, and 
the criticality is explored in the next section. 
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2.3. Analysis of the amplitude ratio m#/A 


In this section, encouraged by the findings in Sec. 12.21 we turn to the 
analysis of the amplitude ratio m H / A, Eq. (P). 

In Fig. [6j we present the scaling plot, (Jnn — ^ArAr)E 1,/ ^-m^/A, for N = 
10,12,..., 22; here, the scaling parameters, J'^ N and u, are the same as those 
of Fig. [H In the XY phase, Jnn — Jff N > 0, the amplitude ratio exhibits 
a plateau for an appreciable range of J NN . Such a feature clearly indicates 
that the amplitude ratio is a universal constant in this domain. 

Upon close inspection, the plateaux in Fig. [6] are curved concavely. The 
shallow bottom locates at J NN = J NN , satisfying dj NN (m H /A)\ j nn =J nn = 0 
for each system size N. We regard the bottom height 

m H /&\j NN =j NN , (7) 

as an indicator for mu/ A. The amplitude ratio (J7J) is plotted for 1 /L 2 [1V(= 
L 2 ) = 10,12,..., 22] in Fig. P The least-squares fit to the data yields an 
estimate m#/A = 2.119(13) in the thermodynamic limit L —>■ oo. As a 
reference, a similar analysis was performed with the abscissa scale replaced 
with 1/L, and we arrived at m#/A = 1.923(17). The discrepancy ~ 0.2 
between these estimates appears to dominate the least-squares-fit error « 
0.02, and the discrepancy may indicate an ambiguity as to the extrapolation 
(systematic error). Regarding it as a possible systematic error, we estimate 
the amplitude ratio as 

m H /A = 2.1(2). (8) 

A comment may be in order, the series of data in Fig. P appear to be 
oscillatory; actually, we observe a slight bump around N(— L 2 ) ~ 16. Such 
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an oscillatory behavior is an artifact of the screw-boundary condition 29], 
rendering an ambiguity as to the extrapolation to L —y oo. The ambigu¬ 
ity appears to be bounded by the above-mentioned error margin, which is 
estimated by performing two independent extrapolation schemes. 


3. Summary and discussions 

The critical behavior of rrifj.G was investigated for the two-dimensional 
quantum S — 1 XY model (J2J) by means of the numerical diagonalization 


method 


29|, 


28]. The interaction parameters were adjusted to Eq. 


m 


order to suppress corrections to scaling 


28]. As a consequence, the data 


(Figs. S| and EJ) collapse into the scaling curves satisfactorily, indicating that 
the data already enter the scaling regime. Thereby, we confirm a universal 
character for the mass-gap ratio (Fig. El), and estimate the amplitude ratio 
as mnl A = 2.1(2). 

As mentioned in Introduction, the amplitude ratio has been estimated 
with the (quantum) Monte Carlo method, m#/A = 2.1(3) [22, 
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24] . as well as the renormalization-group approaches, 2.4 [25], 2.2 26], and 


23] and 3.3(8) 


1.67 [27]. According to the Ginzburg-Landau (mean-field) theory, the am¬ 
plitude ratio should be m H /A = y/2. Clearly, the spectral property reveals 
a notable deviation from that anticipated from the mean-held theory; the 
Ising counterpart was studied in Ref. [32]. In this respect, detailed analyses 


of other spectral properties such as the AC conductivity 23] 133] would be 


desirable. A progress toward this direction is left for the future study. 
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Figure 1: A schematic phase diagram for the two-dimensional 5=1 XY model @ is pre¬ 
sented; here, the parameter space is described by the formula Q. As the nearest-neighbor 
ferromagnetic interaction Jnn increases, a phase transition from the paramagnetic phase 
to the XY phase occurs. The critical behaviors for the Higgs- and Goldstone-mass gaps 
are the main concern. 
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and v = 0.6717 301. 
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( J NN' J Nn)^ 1/V 


Figure 5: The scaling plot, (Jnn — i s presented for the system sizes (+) 

N = 10, (x) 12, (*) 14, (□) 16, (■) 18, (o) 20, and (•) 22. The scaling parameters, J^ N 
and v, are the same as those of Fig. [4j 
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Figure 6: The scaling plot, ( Jnn — J^ N )L X ! V -mu/ A, is presented for the system sizes (+) 
N = 10, (x) 12, (*) 14, (□) 16, (■) 18, (o) 20, and (•) 22. The scaling parameters, J^ N 
and v, are the same as those of Fig. [4] A plateau appears in the XY phase, indicating a 
universal character of the amplitude ratio. 
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Figure 7: The mass-gap amplitude ratio mii /A 0 is plotted for 1/L 2 (N = 10,12,, 22). 
The least-squares fit to the data yields mn /A = 2.119(13) in the thermodynamic limit 
L —> oo. A possible extrapolation error is considered in the text. An oscillatory deviation 
(slight bump around L 2 « 16) is an artifact of the screw-boundary condition 291. 
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